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ABSTRACT 

The Nimbus 7 SBUV measures the same latitude ozone at 
widely different sun angle conditions at the ascent and decent 
part of the orbit during the summer solstice. This situation is 
used in a particular procedure (Ascent/Descent) to obtain the 
relative channel-to-channel calibration error for channels 273 
nm to 306 nm. These estimated errors are combined with 
results from the Pair Justification procedure to correct the 
sun-view diffuser drift in calibration from November 1978 to 
February 1987 for the shorter wavelength channels that mea- 
sure upper stratospheric ozone. Some preliminary re-calibrat- 
ed Nimbus 7 SBUV data in 1989 is compared with the first 
set of SSBUV measurements flown on the Space Shuttle. 


1. INTRODUCTION 

The 1988 ozone trends panel report found the Version 5 
Nimbus 7 SBUV and TOMS total column ozone trends from 
1979 to 1987 to be 3.5 percentage points more negative than 
the Dobson measurements. The Nimbus SBUV profile shows 
a negative 20 percent trend at 48 kilometers. The panel con- 
cluded that the SBUV and TOMS trends were in error due to 
an unaccounted degradation of the instruments’ diffuser 
(Fleig, 1990). Since this finding, a new procedure (designa- 
tion Pair Justification) has been developed and published to 
estimate and correct the diffuser calibration drift for both the 
Nimbus 7 TOMS and SBUV total ozone. TOMS data have 
been reprocessed and released as Version 6 data (Herman, 
1991 ). This paper describes a separate new procedure to esti- 
mate and correct the diffuser calibration drift for the Nimbus 
7 SBUV upper stratospheric ozone (profile) channels. 

The SBUV type instrument measures earth radiance and solar 
irradiance at 12 channels with wavelengths from 255 nm to 
340 nm. Albedos (radiance/irradiance) from 273 nm to 306 
nm are in general used to derive profile ozone. The 255 nm 
channel is not used in Nimbus 7 because of nitric oxide emis- 
sions. The four longest wavelength channels (313, 318, 331 


and 340 nm) are paired to derive total column ozone; A-pair 
(313 and 331 nm), B-pair (318 and 331 nm) and C-pair (331 
and 340 nm). 

The Pair Justification procedure uses a new pair of channels, 
306 nm and 313 nm, designated the D-pair, for assessing total 
ozone. In general, the 306 nm channel is a profile channel. 
However, at equatorial latitudes with low total column ozone 
and small sun angles, 306 nm "sees” the same column ozone 
as the operational pairs (A, B and C). Pair Justification is 
based on a particular comparison of D versus A-pair albedos. 
Analytical leverage for the procedure exists because the D- 
pair is significantly less sensitive to a wavelength dependent 
calibration drift than the A-pair (smaller wavelength separa- 
tion between the channels), plus the D-pair has a greater sen- 
sitivity to ozone. Nimbus 7 TOMS does not have a 306 nm 
channel. Thus, the TOMS' calibration drift is corrected on the 
SBUV D-pair using coincident TOMS nadir-view albedo 
measurements. 

The Pair Justification procedure estimates calibration drift for 
the Nimbus 7 SBUV channels 306 nm to 340 nm. When 
combined with the Ascent/Descent Procedure using a particu- 
lar comparison of ascent versus descent profile measurements 
at or near the summer solstice, calibration drift for channels 
273 nm to 302 nm are estimated and corrected to produce 
Version 6 SBUV ozone data (upper level profile and total 
ozone). 


2. ASCENT/DESCENT PROCEDURE 

During or near summer solstice at the higher latitudes, ozone 
is measured twice in the same orbit, once in the ascent and 
once in descent. Thus, the same latitude is measured twice at 
widely different solar zenith angles. Figure 1 shows the alti- 
tude of normalized single scattered albedo sensitivity to 
ozone for two measurements at the same latitude (72 degrees 
north) with a solar zenith angles of 52 degrees for the ascent 
and 87 degrees for descent. Also shown is the vertical distri- 
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AtrueUsecnt) = ALB[OZRET(A truc (descent))] (4) 

Log[A mea (ascent)] = Log[A true (asccnt)] + E (5) 

Log[A mea (descent)] = Log [A true (dc scent)] + E ( 6 ) 

where A lrue is column matrix of “true" albedos divided by 
albedos for the reference standard profile, A mea is the mea- 
sured albedos divided by albedos for the reference standard 
profile and E is the calibration error. Equation 4. 5 and 6 are 
combine to derive and expression for the calibration error (E). 
A computer code to estimate E uses routines from the Version 
5 Profile Processing Algorithm for the OZRET and ALB cal- 
culations and sets of average ascent and corresponding decent 
measurements (two week, two degree latitude averages). 

Fig. 1. Altitude of normalized single scattered albedo sensitivity to ozone 

for two measurements at 72 degrees latitude during summer solstice plus To test and characterize the computer code formulation, the 

the respective ozone density distribution (normalized). theoretical albedo calculation (ALB) is used to compute sim- 

ulated albedos for a defined ozone profile and simulated cali- 
bution of ozone. In Figure 1 the altitude ot the 273 nm chan- bration error. In other words, the simulated albedos are theo- 

nel (the shortest channel) in ascent is very close to the altitude retical albedos computed from a known ozone profile and cal- 

of the 292 nm channel (a middle channel) in descent. At the ibration error. In instead of the measured albedos, these sim- 

higher solar zenith angles, the altitude where a profile channel ulated albedos are used as input to the computer code formu- 

is sensitive to ozone is displaced upward. Assuming an exact lation. The computed calibration error is compared with the 

overlap of these albedo sensitivities to ozone, a difference in simulated calibration error. Figure 2 shows the results using a 

measured albedo trom a reference standard (Aj^ ea ) versus simulated albedos computed with a calibration error that is a 

calibration error (E) can be described by two equations for the linear function of channel wavelength and ozone profiles 

ascent and descent channels, respectively, as follows. which are uniformly less (same fractional decrease at all alti- 

tudes) than the standard reference. When an ozone profile is 
Log(A m ea273^ = ^273 x ^^true + ^273 ^ the same as the reference standard profile used in the formu- 

lation, there is no difference in the calibration error computed 
Log(A nr)ea 292) = $292 x O^true + ^292 ® by ^ matrix formulation and the simulated calibration error. 

However, as Figure 2 shows, when the profile is different 
where S is the logarithm of albedo sensitivity to changes in from the standard reference, the computed calibration error is 

ozone and OZ^g is the difference in true or actual ozone no t the same as the simulated. On the other hand. Figure 2 

in the region of overlap from the reference standard. With a 
“close" overlap of albedo sensitivity to ozone ascent S 273 IS 
essentially equal to descent S 292 and Equation 1 and 2 
reduces to: 

e 292 ‘ e 273 = L oe( A mea292> ' Lo ?( A mea273) ( 3 ) 

where £ 29 ^ - E 273 * s l ^ e calibration error of 292 nm relative 
to 273 nm. The analytical leverage for computing a relative 
channel-to-channel error from Equation (3) is the result of 
different channels “seeing" the same ozone. As a practical 
matter, the above formulation is not completely independent 
of profile ozone. In order to characterize the altitude of the 
sensitivities, an assumed reference standard ozone profile is 
required. Furthermore, a particular set of ascent and descent 
measurements at a particular latitude band may or may not 
have a set of “very closely" overlapping sensitivities. Thus, a 
matrix formulation is derived which includes an ozone profile 
retrieval (OZRET) and theoretical albedo calculation (ALB) 
and the following overall relationships: 



Fig. 2. Computed calibration errors using simulated albedos from a sim 
ulated calibration error which is a linear function of wavelength and from 
ozone profiles that are uniformly less than the reference standard profile. 
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shows that the relative channel-to-channel calibration error is 
accurately determined. Figure 3 shows the results for a simu- 
lated calibration error that is a non-linear function of channel 
wavelength and “bump” changes in profile ozone. A “bump” 
has minus 20 percent change in ozone in the designated 
Umkehr Layer and minus 10 percent change in the layer 
above and the layer below. As with the uniform differences 
in ozone in Figure 2, the “bump” ozone results in large differ- 
ences in computed versus simulated calibration error but the 
relative channel-to-channel calibration error has useful accu- 
racy. 

3. CALIBRATION DRIFT FOR PROFILE CHANNELS 

Figure 4 shows how the relative channel-to-channel calibra- 
tion error computed by the formulation can be combined with 
other procedures to develop a “total calibration error” for all 
channels. Total calibration error includes both a “Drift” and 
an “Initial” calibration error. As discussed above, to-date 
only the drift in calibration for the longer total ozone wave- 
lengths has been determined by Pair Justification. Thus, the 
drift in the computed relative channel-to-channel calibration 
error for channels 273 to 306 nm is joined to the results of 
this Pair Justification to estimate the calibration drift for all 
the channels. A calibration correction for this estimated drift 
is used to produce Version 6 Nimbus 7 SBUV ozone profiles. 


4. NIMBUS 7 SBUV CALIBRATION DRIFT RESULTS 

For Nimbus 7 SBUV, only the northern hemisphere summer 
solstice data is used. In the southern hemisphere, there 
appears to be a instrument temperature anomaly as the satel- 
lite comes out of darkness causing some differences between 
the ascent versus the descent orbit measurements. Work is 
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Fig. 3. Computed calibration errors using simulated albedos from a sim- 
ulated calibration error which is a non-linear function of wavelength and 
from ozone profiles that have "bump ' reductions in ozone from the refer- 
ence standard profile. 


Fig. 4. Overall procedural steps to derive an estimated total calibration 
error. 

planned to further investigate this issue. In the northern 
hemisphere, individual scan albedo data is averaged for two 
degree latitude bands during two week periods. Figure 5 
shows the computed Version 5 relative channel-to-channel 
calibration error from six sets of these averages during the 
1986 northern summer solstice. As with the simulation, there 
are large differences between the different averages but a con- 
sistent relative channel-to-channel error. The drift in calibra- 
tion error is computed by subtracting the result of each year 
from the average of years 1980 and 1981. Figure 6 shows the 
computed drift in 1986 from the average of the same six data 
set in Figure 5. Also shown is the drift for the longer wave- 
length channels from Pair Justification. The relative channel- 
to-channel calibration error from 273 nm to 306 nm is con- 
nected to the Pair Justification results. Figure 7 shows the 
resulting estimated drift for all channels for the month of June 
from 1985to 1989. 



Fig. 5. Computed relative channel-io-clKinnel calibration error from 
two-week, two-degree-laiiludc average albedos. 
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Tig 6. Drill in relative channel- (o-channcl calibration error versus drift 
in longer wavelengths from pair justification. 


Starting in March 1987, the Nimbus 7 SBUV measurements 
developed a large increase in noise due to non-synchroniza- 
tion of the chopper wheel. The evidence suggest there are 
greater inaccuracies in estimating the calibration error for this 
data. A procedure to correct the non-synchronization is 
presently under development. Thus, the data shown in Figure 
7 for 1987 through 1989 is preliminary and has larger uncer- 
tainties than the pre-March 1987 data. 

For the pre-March 1987 data, the largest uncertainty in the 
estimated calibration drift for 273 nm is the uncertainty in 
306 nm from Pair Justification. The second- largest uncertain- 
ty is the uncertainty in solar zenith angle determination. At 
mid-latitudes, the combined uncertainties, gives an estimated 
uncertainty in ozone derived from the 273 nm channel of ±6.4 
percent ozone (95 percent confidence). For ozone derived 
from 306 nm channel, the uncertainty is slightly less at ±5.0 
percent ozone. 



The SSBUV instrument flown on the Space Shuttle is 
designed to validate the calibration of SBUV type instru- 
ments. Its first Bight was in October 1989. Figure 9 shows 
the average differences in SSBUV versus Nimbus 7 Version 5 
(Uncorrected) SBUV coincident albedos. Also shown is a 
preliminary estimated calibration drift from June 1979 to 
October 1989. However, the comparison of SSBUV and the 
estimated drift are not fully complimentary. As shown in 
Figure 4, the total calibration error also needs an estimate of 
the “Initial Error”. In the case of Nimbus 7, the first summer 
solstice is eight months after launch. Thus, an estimate of the 
“Initial Error” using June 1979 data includes both an estimat- 
ed “at launch” calibration error and the drift from launch to 
June 1979. The procedure to estimate “Initial Error” is 
presently under development. 



Fig. 8. Comparison of total calibration error in October 1089 estimated 
by SSBUV versus estimated calibration drift by Ascent/Descent Procedure. 

5. CONCLUSION 

The drift in Calibration error for Version 5 Nimbus 7 SBUV 
has been corrected for data from November 1979 to February 
1987. The uncertainty in mid-latitude profile ozone derived 
by the shortest wavelength channel is ±6.4 percent (95 per- 
cent confidence). Preliminary data for total calibration error 
for October 1989 shows a relatively close comparison with 
the first SSBUV Bight 
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Fig 7 Estimated calibration dnft for all channels used for Version 6 
calibration correction for the month of June, years 1085 through 1080; 
ftelimituiry result for 1087 through 1080. 
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